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SUMMARY 

A modular  model  of  the  sudden-expansion  combustor  has  been 
developed,  combining  a parabolic,  directed-f low  computation 
and  a well-stirred  reactor  model  representing  the  recircula- 
tion region(s).  The  parabolic  module  incorporates  a two-equation 
turbulent  kinetic  energy  turbulence  model,  and  both  modules  use 
the  quasiglobal  formulation  for  rapid  computation  of  finite-rate 
chemical  reactions.  Coupling  between  the  modules  is  through  a 
shear  layer  representation  which  provies  the  parabolic-flow 
boundary  conditions  and  the  stirred-reactor  feed  rates.  Compared 
to  unified  elliptic  models  of  the  dump  combustor,  the  major 
advantage  of  this  approach  is  the  greatly  increased  detail  possible 
in  the  combustor  analysis.  This  detail  in  turn  allows  the  develop- 
ment of  chemical  kinetic  and  flowfield  models  applicable  to  the 
refinement  of  elliptic  formulations. 
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ABSTRACT 

This  report  describes  the  development  of  a modular  model 
for  the  prediction  of  the  performance  of  sudden  expansion  burners 
as  a function  of  the  controllable  parameters  relevant  to 
combustor  design.  The  model  is  based  upon  a concept  in  which  the 

recirculation  zone,  treated  as  a stirred  reactor,  is  coupled 
to  a parabolic  boundary  layer  formulation  for  the  flow  outside 
the  recirculation  zone.  Hydrocarbon  oxidation  kinetics  and 
turbulent  kinetic  energy  models  are  employed  in  the  model  develop- 
ment. In  addition  to  the  parabolic-flow  and  stirred  reactor 
elements,  a module  representing  the  fuel  injection  process  has 
been  developed.  Results  of  the  application  of  the  modular  model 
to  the  analysis  of  cold-flow  and  reacting-f low  dump  combustor 
experimental  data  are  described. 
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NOMENCLATURE 

MEANING 

Coefficient  in  shear  layer  width  expression 
(growth  rate) 

Coefficient  in  shear  layer  width  expression 
(initial  thickness) 

Coefficients  in  two-equation  turbulence 
model 

Skin  friction  coefficient 

Prandtl  model  eddy  viscosity  coefficient 

Combustor  diameter 

Dump  height  = (D3~D2)/2 

Turbulent  kinetic  energy 

Shear  layer  width 

Turbulent  kinetic  energy  dissipation  length 
scale 

Width  scale  in  Prandtl  eddy  viscosity  model 
Mach  number 
Total  pressure 
Static  Pressure 

Prandtl  number  for  turbulent  kinetic  energy 

Radial  coordinate;  with  subscript  radius 

Distance  from  wall  to  first  parabolic 
flow  grid  point 

Total  temperature 
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MEANING 

Axial  mean  velocity  component 
Mass-average  axial  mean  velocity 
Radial  mean  velocity  component 
Axial  coordinate 

Turbulent  kinetic  energy  dissipation  rate 

Turbulent  dynamic  viscosity 

Density 

Turbulent  shear  stress 
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Centerline 

Dividing  streamline  (boundary  between 
directed  flow  and  stirred  reactor  in 
modular  model) 

Grid  point  in  parabolic  flow  just  inside 
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Recirculation  zone  value 
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1.  INTRODUCTION 

Current  and  future  airbreathing  propulsion  system  concepts 
are  tending  toward  lower  weight  and  smaller  volume  designs. 

An  example  is  the  integral  rocket/ramje  engine  which  uses  a 
common  combustor  for  rocket  boost  and  sustained  ramjet  flight. 
Unfortunately,  the  current  design  trends  create  stringent  re- 
quirements for  higher  performance  from  smaller  volumes,  but  the 
combustion  processes  are  not  well  enough  understood  to  optimize 
combustor  designs.  Thus,  combustor  development  has  been  carried 
out  largely  by  employing  costly  cut-and-try  methods. 

The  basic  difficulty  in  understanding  combustor  behavior 
is  the  coupling  between  turbulent  mixing  and  reaction  kinetics. 
This  problem  is  not  unique  to  combustor  flows,  but  is  encountered 
in  other  systems  where  flow,  mixing,  and  reaction  times  are 
all  of  the  same  order  of  magnitude,  for  example,  external  burn- 
ing devices  for  thrust  and/or  control  and  gas  dynamic  and 
chemical  laser  systems.  In  addition,  many  of  the  reacting 
flows  of  general  interest  involve  condensed  phase  species  and 
large  scale  recirculation  zones  such  as  are  found  in  the  "dump" 
burner  concept.  Problems  in  ignition,  flame  stabilisation, 
combustion  efficiency,  and  flow'field  and  combustion  instabilities 
also  occur  in  the  dump  burner  as  well  as  in  many  other  combustion 
systems  of  interest,  and  these  problems  need  to  be  better  under- 
stood to  support  engineering  judgement. 

Consequently,  continued  research  is  required  in  both  sub- 
sonic and  supersonic  flows  in  order  to  gain  the  necessary  under- 
standing of  the  controlling  mechanisms.  Only  in  this  way  can 
effective  design  principles  be  developed  to  provide  optimum 
propulsion  system  performance  over  wide  operating  ranges. 

The  objective  of  the  work  summarized  in  this  report  is  to 
eliminate  the  deficiencies  just  described  and  provide  a sound 
basis  for  effective  combustor  design  criteria,  through  the 
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development  of  analytical  models  for  the  sudden-expansion  comb- 
ustor flowfield.  The  major  emphasis  in  this  program  is  on  the 
development  of  a "modular"  model  for  the  sudden  expansion  comb- 
ustor, in  which  proved  computational  elements,  including  a 
parabolic  flowfield  analysis  procedure  and  a perfectly-stirred 
reactor  formulation,  are  coupled  together  to  model  the  overall 
flowfield.  In  this  report  the  features  of  the  modular  approach 
are  reported  and  progress  in  the  development  and  application 
of  the  model  to  cold-flow  and  reacting-flow  combustor  data  is 
described.  This  work  has  involved,  and  will  continue  to  in- 
volve, close  coordination  with  the  experimental  work  in  progress 
at  the  Air  Force  Aero  Propulsion  Laboratories  at  Wright-Patterson 
Air  Force  Base,  as  well  as  other  related  programs. 
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2.  MODULAR  MODEL  CONCEPT 

The  basic  characteristics  of  the  sudden  expansion  burner 
are  depicted  in  Figure  1.  In  essence,  the  key  feature  of  the 
flow  field  is  the  recirculation  region  just  downstream  of  the 
expansion,  separated  from  the  directed  flow  by  the  dividing 
streamline.  Because  of  this  recirculation  region,  in  which 
there  is  no  characteristic  flow  direction,  the  flow  field  in 
the  sudden  expansion  burner  is  in  general  described  by  ellip- 
tic equations;  i.e.,  the  boundary  layer  approximations  cannot 

be  applied  to  the  whole  flow  and  the  elliptic  Navier-Stokes  l 

equations  must  be  used.  However,  it  is  possible,  at  least 
conceptually,  to  divide  the  flow’  field  into  two  regions.  The 
first,  the  directed  flow  region,  consists  of  that  part  of  the 
flow  between  the  dividing  streamline  and  the  burner  centerline. 

The  second  region  is  the  recirculating  flow  region,  and  the  two 
regions  are  coupled  through  the  conditions  at  the  dividing 
streamline . 

The  concept  of  the  modular  approach  to  the  computation  of 
the  flow  in  a sudden  expansion  burner  involves  the  application 
of  available  techniques  to  the  separate  regions  of  the  flow- 
field.  Thus,  the  directed  flow  region  is  treated  using  a para- 
bolic finite  difference  solution  of  the  boundary  layer  form  of 
the  equations  of  motion.  The  recirculation  regions  are  treated 
as  perfectly  stirred  reactors,  with  the  reactor  feed  rate 
defined  by  the  species  diffusion  flux  across  the  dividing  stream- 
line. The  shape  of  the  dividing  streamline  must  be  specified 

I 

a priori;  the  resulting  wall  static  pressure  distribution  may 
then  be  compared  with  experiment  and  another  dividing  streamline 
shape  specified,  if  necessary.  Since  the  flow  field  is  in  general 
elliptic,  the  modular  approach  approximates  the  true  flow  field, 
and  the  specification  of  the  species  fluxes  across  the  dividing 
streamline  must  be  carried  out  iteratively.  However,  despite 


i 


3 


Figure  1.  Schematic  of  Sudden  Expansion  (Dump) 
Burner 
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its  limitations,  this  method  does  allow  the  application  of 
well-tested  and  proved  flow  models  to  the  sudden  expansion 
burner  geometry,  and  is  capable  of  yielding  parametric  infor- 
mation of  interest  in  sudden  expansion  burner  design. 
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3. MODULAR  MODEL  FORMULATION 


The  basic  elements  of  the  current  modular  model  formulation 
are  the  parabolic  finite-difference  computational  technique 
developed  by  Boccio,  Weilerstein  and  Edelraan  [l],  modified  to 
incorporate  the  turbulent  kinetic  energy  model  developed  by 
Harsha  [2],  used  for  the  directed-f low  portion  of  the  analysis, 
and  the  stirred  reactor  computation  developed  by  Edelman  and 
Weilerstein  [3] . Both  of  these  elements  make  use  of  the  quasi- 
global  model  developed  by  Edelman  and  Fortune  [4]  for  rapid 
computation  of  finite-rate  hydrocarbon  air  kinetics.  These 
elements,  or  modules,  are  coupled  together  through  a simplified 
representation  of  the  turbulent  shear  layer  which  exists  between 
the  directed  flow  and  the  recirculation  region.  In  the  modular 
approach,  the  shear  layer  representation  is  used  to  define  the 
gradients  in  velocity,  species,  and  enthalpy  between  the  two 
regions  of  the  flow,  thus  providing  both  the  boundary  conditions 
on  the  directed  flow  and  the  stirred  reactor  feed  rates.  Details 
of  the  mathematical  formulations  used  for  the  flowfield  region 
are  discussed  in  this  section;  details  of  the  numerical  analysis 
procedure  may  be  found  in  Ref.  [5]. 

A detailed  description  of  the  mathematical  formulation  of 
the  modular  model  has  been  given  in  a previous  AFOSR  interim 
report,  Ref.  6 , in  which  results  of  the  application  of  the 
model  to  the  prediction  of  both  cold  flow  and  reacting  flow  in 
a sudden-expansion  combustor  were  reported.  One  observation  from 
these  previous  results  was  that  the  algebraic  length  scale 
formulation  used  in  the  one-equation  turbulent  kinetic  energy 
model  developed  by  Harsha  [2]  was  inadequate  for  sudden  expansion 
combustor  geometries.  The  reason  for  this  inadequacy  is  that 
there  are  two  major  characteristic  length  scales  that  can  be 
identified  in  the  sudden-expansion  combustor  flowfield,  one 
applying  to  the  main  flow  (directed  flow)  and  one  to  the  shear 
layer  region.  The  radial  variation  in  length  scale  across  the 


i 
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combustor  that  develops  from  the  existence  of  these  two  characteristic 
length  scales  cannot  be  easily  modeled  using  an  algebraic  length 
scale  formulation. 

To  provide  a means  of  accounting  for  radial  length  scale 
variation  in  the  sudden-expansion  combustor  geometry,  the  two- 
equation  turbulence  model  of  Launder,  et  al.  [ 7]  was  incorporated 
into  the  directed  flow  portion  of  the  modular  model.  This  turbulence 
model  involves  the  solution  of  transport  equations  for  the  turbulent 
kinetic  energy  and  the  turbulence  energy  dissipation  rate: 


Turbulent  Kinetic  Energy 

8k  3k  _1  3,yTr  3k.  . ,3u2 

Pu  3x  3 j*  j.  g"j£pr  3r ^ yT^3r  ~ P^ 

k 

Turbulence  Energy  Dissipation 


(1) 


3e 


. 3 c 13  yTr  3e . . _ e ,3u? 

pu  3x  + pv  ITr  “ r ITr  (PFT  37)  + Cel  k yT(3?) 


Ce2pl 


(2) 


where 


yT  = Cyp  k2/e  (3) 

and  in  which 

Cy  = 0.09,  C£l  = 1.40,  Ce2  = 1.95 

The  turbulence  energy  dissipation  rate,  e,  can  be  related  to  a 
turbulent  length  scale,  through  the  Kolmogorov  hypothesis 

e <*  k3/2/£k 


so  that  the  transport  equation  for  turbulence  energy  dissipation 
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can  also  be  regarded  as  a description  of  the  spatial  variation 
of  the  turbulent  length  scale. 


Boundary  conditions  for  the  turbulent  kinetic  energy,  k, 
and  the  turbulent  kinetic  energy  dissipation  rate,  e,  are  established 
in  the  same  manner  as  described  in  Ref.  [6]  for  the  mean  total 
energy  and  species  equations.  Thus,  values  are  assumed  for  the 
turbulent  kinetic  energy  and  dissipation  length  scale  for  the 
stirred  reactor  region,  k^  and  respectively,  and  the  gradient 

in  k and  e along  the  dividing  streamline  established  from  the 
relations 


3k  _ ki'  krZ 
3r  i 


and 

3e  _ ei  ~ eRZ 
3r  l 


where  the  subscript  i represents  the  first  grid  point  inside  the 
directed  flow  region  module  and  £= ax+b  is  the  assumed  shear 
layer  width. 


It  should  be  pointed  out  that  this  approach  is  substantially 
different  from  that  described  in  Ref.  [6]  for  the  modular  model 
using  a one-equation  turbulent  kinetic  energy  formulation.  For 
the  one-equation  model,  the  value  of  the  turbulent  kinetic  energy 
along  the  dividing  streamline  was  related  directly  to  the  dividing 
streamline  shear  stress,  obtained  from  the  assumed  dividing  stream- 
line skin  friction  coefficient.  While  the  approach  used  in  the 
present  formulation  introduces  two  additional  parameters,  k^z 
and  , it  has  been  found,  as  will  be  shown  in  the  next  section, 
that  this  method  results  in  a considerably  reduced  sensitivity 
of  the  prediction  to  the  values  specified  for  the  coefficients 
a and  b in  the  shear  layer  width  expression.  Furthermore,  this 
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approach  allows  the  investigation  of  the  sensitivity  of  the 
overall  flowfield  prediction  to  the  "intensity"  of  the  mixing 
in  che  stirred  reactor  region,  as  measured  by  and 

lor  combustor  regions  in  which  the  parabolic  directed  flow 
computation  is  carried  to  the  wall,  boundary  conditions  on  the 
turbulent  kinetic  energy  and  its  dissipation  rate  are  established 
in  a manner  similar  to  other  applications  of  the  two-equation 
model  to  ducted  flow.  Thus,  the  wall  value  of  kinetic  energy  is 
directly  proportional  to  the  wall  shear  stress,  obtained  from 
the  specification  of  a wall  skin  friction  coefficient,  and  the 
value  of  the  dissipation  rate  at  the  wall  is  obtained  from  the 
Kolmogorov  hypothesis  using  the  wall  turbulent  kinetic  energy 
and  a local  length  scale  related  to  the  law-of-the-wall  for  a 
turbulent  boundary  layer.  It  should  be  noted  that  within  the 
context  of  the  modular  model  for  sudden  expansion  combustors  the 
effects  of  the  specification  of  the  flowfield  boundary  conditions 
at  the  combustor  wall  are  negligible:  it  is  the  specification 
of  boundary  conditions  along  the  dividing  streamline  that  dominates 
the  solution  to  the  problem. 
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4.  APPLICATION  OF  THE  MODULAR  MODEL  TO  PREDICTION 
OF  THE  PERFORMANCE  OF  A PROTOTYPE  SUDDEN-EXPANSION  COMBUSTOR 

A comprehensive  study  of  the  performance  of  center  dump 
combustors,  with  and  without  f lameholders , has  been  reported  by 
Craig,  et  al.  (Ref.  [8]  ).  This  paper  includes  results  for  both 
premixed  and  wall  fuel  injection  configurations.  In  the  wall 
fuel  injection  configuration,  fuel  was  injected  through  8 equally 
spaced,  0.14  cm  dia.  orifices  located  12.07  cm  upstream  of  the 
dump  plane,  while  in  the  premixed  configuration  fuel  was  injected 
through  spraybars  located  1.5  m upstream  of  the  dump  plane.  For 
all  of  the  results  reported  in  Ref.  [8]  , the  fuel  was  JP-4 . 

The  basic  configuration  studied  in  Ref.  8 involved  a 
combustor  L/D  of  3,  a dump  area  ratio  A^/Ag  of  2.25,  and  an  exit 
nozzle  area  ratio  Ag/A^  = 0.40.  Two  sets  of  inlet  airflow  rates 
and  inlet  total  temperatures  were  used  throughout,  with  nominal 
values  of  1.45  kg/sec  and  0.73  kg/sec,  and  555°K  and  722°K, 
respectively.  A major  contribution  of  this  work  is  the  measurement 
of  wall  static  pressure  distributions  in  both  cold  flow  and  reacting 
flow.  This  report  represents  one  of  the  few  investigations  in 
which  these  data,  crucial  to  the  comparison  of  experimental  data 
with  analytical  model  results,  have  been  reported.  Supplementing 
the  results  reported  in  Ref.  [8]  , additional  unreported  static 
pressure  distributions  have  been  made  available  to  the  authors 
by  Dr.  Craig. 
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4.1  EFFECTS  OF  MODEL  COEFFICIENTS  IN  COLD  FLOW 

As  is  obvious  from  the  description  of  the  modular  model 
in  Section  3 and  in  Ref.  [6]  , there  are  a number  of  parameters 
in  the  model  which  must  be  specified  in  order  to  compute  a given 
flowfield.  These  include  geometric  parameters  such  as  the  size 
and  shape  of  the  recirculation  zone,  and  the  initial  thickness 
and  spread  rate  of  the  shear  layer  along  the  dividing  streamline 
between  the  recirculation  zone  and  the  directed  flow,  and  flow- 
field  parameters  including  the  effective  turbulent  intensity  and 
scale  within  the  recirculation  zone  and  the  effective  turbulent 
shear  stress  along  the  dividing  streamline.  Thus  it  is  of  considerable 
importance  to  determine  the  sensitivity  of  the  model  predictions 
to  these  parameters,  as  well  as  the  sensitivity  of  the  model 
predictions  to  inlet  condition  specification  and  to  the  coefficients 
used  in  the  turbulence  model  itself. 

Several  of  the  model  coefficients  are  related  to  each  other, 
although  it  is  not  always  possible  to  establish  the  form  of  the 
relationship  in  an  a priori,  general  fashion.  Thus,  the  dividing 
streamline  turbulent  shear  stress  level  can  be  expected  on  physical 
grounds  to  be  related  to  the  shear  layer  growth  rate  and  the 
effective  turbulent  intensity  and  length  scale  level  in  the 
recirculation  zone,  for  example.  However,  for  purposes  of  this 
evaluation  of  the  effects  of  model  coefficients,  each  of  the 
parameters  was  varied  independently.  Furthermore,  it  is  of 
considerable  interest  from  the  standpoint  of  model  evaluation 
to  determine  the  accuracy  with  which  a reacting  dump  combustor 
flowfield  can  be  predicted  using  model  coefficients  established 
by  comparison  with  a corresponding  cold  flow  experiment,  both 
to  evaluate  the  effects  of  reaction  on  the  overall  flowfield  and 
because  more  detailed  data  are  generally  available  for  nonreacting 
flowfields.  The  data  reported  by  Craig,  et  al.,  Ref .[ 8] , includes 
a broad  range  of  experimental  conditions  for  geometrically  similar 
combustors,  in  6oth  cold  flow  and  hot  flow,  and  includes  wall 
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static  presure  distributions,  which  earlier  work  (Ref.  6 ) has 

indicated  can  be  used  to  assess  the  effects  of  parameter  variations 
on  the  model  prediction. 

For  the  cold  flow  calculations,  the  cold  flow  data  obtained 
by  Craig  et  al . (Ref.  [8])  in  an  area  ratio  2.25  center-dump  combustor 
geometry  at  an  inlet  total  temperature  of  552°K  and  an  air  mass 
flow  of  1.45  kg/sec  were  used  for  comparison.  The  initial  conditions 
for  the  calculations  were  thus: 

P-  = 2045.3  lb/ft2  (0.97  atm) 

2 

T0  = 930°R  (516. °K) 

S2 

Ug  = 895  ft/sec  (273  m/sec) 

M2  = 0.603 

The  initial  turbulent  kinetic  energy  level  was  taken  to  be 
K2  = 0.03(U2)2 

and  the  initial  turbulent  kinetic  energy  dissipation  rate  was 
obtained  from 

e2  = 1.69(k2)3/2/r2 

An  adiabatic  combustor  wall  was  assumed  in  all  calculations  and 
the  wall  skin  friction  coefficient  (applied  in  the  regions  up- 
stream and  downstream  of  the  recirculation  zone)  was  taken  to 
be  0.003. 

The  values  of  the  parameters  in  the  modular  model  can  in 
most  cases  be  bounded:  for  example,  the  coefficient  b in  the 
shear  layer  width  expression  must  lie  in  the  range  0<b<r2,  while 
similarly  a must  be  bounded  such  that  at  the  end  of  the  recirculation 

zone  0<ax+b<ro-  The  recirculation  zone  turbulent  kinetic  energy 

**  2 
can  be  assumed  to  be  of  the  order  of  magnitude  of  (U2)  , while 
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the  length  scale  which  defines  the  recirculation  zone  turbulent 
kinetic  energy  dissipation  rate  should  be  of  the  order  of  the 
dump  step  height.  Thus  for  the  configuration  used  by  Craig 
et  al.  [8]  we  have 


0<a<0. 32 
0<b<0. 167 

kRz=  <U2>2 

XRz  * h 

An  initial  set  of  calculations  of  the  Craig  et  al.  cold  flow 

data  was  carried  out  using  CpRz>  the  dividing  streamline  "skin 

friction"  coefficient,  as  a parameter,  with  a = .02,  b = 0.07, 

2 

k^  = (U0)  , and  1R  = h/2.  Results  of  these  computations  are 


from  0.12  to  0.18.  As 


Rz  v~2'  ’ Rz 

shown  in  Fig  2a,  for  a range  in  Cppz 

would  be  expected,  an  increase  in  CpRz  results  in  a greater 

velocity  retardation  effect  on  the  directed  flow,  as  is  indicated 

by  the  increase  in  value  of  the  parameter  (ucl~u)/u  shown  in 

Fig.  2b.  Here  u represents  the  mass-average  velocity  in  the 

directed  flow:  for  uniform  flow  (u  ,-u)/u=0,  while  for  a para- 

ci  o . _ 

bolic  velocity  profile,  u/uc^=  l-(r/rw)  , (uc-^-u)/u  = 0.5. 

Although  the  CpRz  parameter  cannot  be  directly  related  a 
priori  to  the  turbulent  viscosity  level  along  the  dividing 
streamline,  since  the  viscosity  is  determined  through  the  computed 
divinding  streamline  turbulent  kinetic  energy  and  dissipation 
rates,  it  is  of  interest  to  compare  the  computed  turbulent  viscosity 
levels  to  that  which  would  be  obtained  using  an  eddy  viscosity 
formulation,  for  example  the  Prandtl  model 


P^PR^  PR 

Using  the  computed  values  of  uT  and  Au,  for  the  CpRz=0.18  case, 
and  taking  lp^  to  be  the  shear  layer  width  scale  1,  results  in 
a range  of  values  of  the  coefficient  Cp  of  0.37<  Cpp<.042  along 
the  length  of  the  shear  layer.  These  are  within  a factor  of  two 
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b)  Effect  of  Recirculation  Zone  "Skin  Friction" 
Coefficient  on  Velocity  Profile  Development 


FIG.  2.  PARAMETER  VARIATION  EFFECTS  ON  MODULAR  MODEL  COLD 
FLOW  PREDICTION,  DATA  FROM  CRAIG  8 , L/D=3 , 

A3/A2=2-25.  A*/A3=0.40,  TqI=  1000°R,  M^O.603 
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of  the  values  of  CpR  commonly  used  for  turbulent  shear  layer 
calculation  using  the  Prandtl  eddy  viscosity  model. 

From  thse  results,  it  is  evident  that  the  best  prediction 
of  the  cold  flow  data  of  Craig  et  al.  [8]  is  achieved  using 
CpRz=0.18.  The  next  parameter  to  be  varied  is  the  initial  r'v>ear 
layer  thickness,  b,  and  the  results  of  variation  of  this  parameter 

over  a range  .05<b<.09  are  shown  on  Fig.  2c.  Here  it  is  evident 

that  the  major  effect  of  the  initial  shear  layer  thickness  is 
on  the  pressure  distribution  in  the  early  portion  of  the  flow, 
just  downstream  of  the  dump  plane  (x/D=0).  Decreasing  the  initial 
shear  layer  thickness  decreases  the  initial  pressure  drop  observed 
both  in  the  data  and  in  the  predictions,  at  a given  value  of  CpRz- 

This  is  consistent  with  a decrease  in  turbulent  viscosity  level 

along  the  dividing  streamline,  which  does  in  tact  occur.  However, 
the  magnitude  of  the  decrease  is  considerably  less  than  would 
be  predicted  by  a simple  eddy  viscosity  model  with  change  in 
length  scale,  because  of  the  use  of  a two-equation  turbulence 
model  to  obtain  the  viscosity.  As  the  shear  layer  thickness 
is  decreased,  the  transport  of  both  turbulence  energy  and  its 
dissipation  rate  from  the  recirculation  zone  into  the  directed 
flow  is  increased  proportionally.  The  net  effect,  however, 
is  to  decrease  the  turbulent  kinetic  energy  along  the  dividing 
streamline,  because  of  an  increase  in  enery  dissipation  at  this 
point  in  the  flow.  This  results  in  a decrease  in  turbulent  viscosity, 
but  a decrease  smaller  than  that  which  would  be  predicted  by 
a simple  viscosity  model.  Consequently,  the  overall  effect  of 
a +30%  change  in  initial  shear  layer  thickness  on  the  predicted 
pressure  gradient  is  relatively  small,  as  shown  in  Fig.  2c. 

Changes  in  the  shear  layer  growth  rate  constant,  a,  produce 
essentially  no  change  in  the  predicted  static  pressure  distribution, 
as  can  be  seen  from  Fig.  2d,  for  the  range  .01<a<03.  Since  the 
shear  layer  length  to  initial  thickness  ratio  for  this  configuration 
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c)  Effect  of  Shear  Layer  Initial  Thickness  on 
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d)  Effect  of  Shear  Layer  Growth  Rate  on 
Wall  Static  Pressure  Distribution 
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(with  b=0.05)  is  15.3,  the  width  of  the  shear  layer  increases 
46%  for  a=0.03  and  only  15%  for  a=0.01,  so  that  the  initial 
thickness  dominates  the  transport  rate  across  the  shear  layer 
over  most  of  its  length.  This  observation,  and  the  relative 
insensitivity  of  the  model  to  shear  layer  thickness,  with  the 
two-equation  turbulence  formulation,  provides  an  explanation 
for  the  insensitivity  observed  in  Fig  2d. 


The  remaining  parameters  involved  in  the  model  formulation 


are  the  recirculation  zone  kinetic  energy, 


k and  the  effective 
rz 


dissipation  length  scale  in  the  recirculation  zone,  1 


rz 


For 


the  computations  so  far  described,  k was  set  equal  to  the 

square  of  the  inlet  velocity,  and  1 was  taken  to  be  h/2, 

where  h is  the  dump  step  height.  Parametric  calculations 

2 2 

with  k ^ as  the  parameter  were  carried  out  for  0.5uo<k  <1.5uo* 

rz  z rz  z 

Note  that  since  the  recirculation  zone  dissipation  rate  e =1.69 

3/2  rz 

(krz)  /lrzi  the  value  of  erz  for  each  of  these  cases  was 

adjusted  accordingly. 

2e . 


Results  of  these  parametric  calculations 
are  shown  on  Fig.  2e.  As  can  be  seen  from  Fig.  2e  the  results 
of  this  parametric  evaluation  show  a greater  initial  pressure 
drop  occurring  as  k decreases.  This  is  indicative  of  an  increasing 
dividing  streamline  turbulent  viscosity  coefficient  as  the 
recirculation  zone  kinetic  energy  level  decreases.  Similarly, 
the  effects  of  variation  of  the  assumed  recirculation  zone 


as  1 is  decreased, 


length  scale,  1 are  shown  in  Fig.  2 f : 
r z 

for  given  k , e increases,  so  that  the  viscosity  level  along 
rz  rz 

the  dividing  streamline  decreases  and  the  initial  pressure  drop 

decreases.  It  should  be  noted  with  regard  to  these  comparisons 

that  the  parametric  variations  involved  in  the  results  decribed 

are  somewhat  artificial,  since  in  both  cases  shown  in  Figs.  2e 

and  2f  the  effective  dividing  streamline  shear  stress  was  held 

constant.  Clearly  the  specification  of  dividing  streamline 

shear  stress  must  affect  the  turbulent  intensity  level  within 

the  recirculation  zone,  since  physically  it  is  the  conversion 
*A  value  of  1.5  u|  is  physically  unrealistic  but  has  been  used 
to  determine  the  sensitivity  of  the  model  to  large  parameter 
variations 


17 


AXIAL  DISTANCE /COMBUSTOR  DIAMETER,  X/D 


e)  Effect  of  Recirculation  Zone  Turbulent 

Kinetic  Energy  Level  on  Wall  Static  Pressure 
Distribution 
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Static  Pressure  Distribution 
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of  directed  flow  mean  kinetic  energy  to  recirculation  zone 
turbulent  kinetic  energy  through  the  shear  layer  region  that 
drives  the  mixing  in  the  recirculation  zone.  This  coupling  has 
not  been  accounted  for  in  the  present  computations,  but  will 
be  modeled  in  further  development  of  the  modular  formulation. 

As  a result  of  these  parametric  calculations,  the  coefficients 

chosen  for  the  reacting  flow  computations  were  a=0.03,  b=0.05, 

2 

CFrz=0-18,  krz=(u2^  and  lrz=h/2.  These  cold  flow  computations, 
besides  establishing  a preliminary  set  of  coefficients  for  the 
reacting  flow  calculations,  also  show  that  the  use  of  the  two- 
equation  turbulence  model  reduces  the  sensitivity  of  the  combustor 
model  to  the  values  chosen  for  the  shear  layer  coefficients  a 
and  b.  The  wall  static  pressure  predictions  are  sensitive  to 
the  values  chosen  for  CFrZ,  k^z , and  1RZ,  with  the  apparent 
sensitivity  involving  both  the  dividing  streamline  viscosity 
implied  by  these  parameters  and  the  relative  transport  rates  of 
turbulent  kinetic  energy  and  its  dissipation  rate. 

4.2  EFFECT  OF  MODEL  COEFFICIENTS  IN  PREMIXED  REACTING  FLOW 

With  the  appropriate  model  coefficients  determined  from  the 
cold  flow  calculations  just  outlined,  predictions  were  made  for 
comparison  with  the  premixed  dump  combustor  data  reported  by 
Craig  et  al.  Ref. [8]  . The  predictions  were  made  for  a fuel- 
air  ratio  of  0.053.  In  the  experiments,  the  fuel  used  was  JP-4 ; 
for  the  computations  propane  was  used  to  represent  the  fuel. 

A simple  one-step  global  finite  rate  chemistry  model  was  in- 
corporated in  these  computations. 

Initial  conditions  included  an  inlet  total  temperature  of 
554°K  (997°R)  and  an  inlet  static  pressure  of  1.83  atm  (3881 
lb/ft^).  With  a mass  flow  rate  of  1.57  kg/sec  (3.45  lb/sec), 
the  inlet  velocity  and  static  temperature  were  159  m/sec  (521 
ft/sec)  and  543°K  (977°R),  respectively;  the  inlet  Mach  number 
was  0.351.  Initial  turbulent  kinetic  energy  and  dissipation 
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rate  values  were  established  in  the  same  manner  as  for  the  cold 
flow  calculations,  and  the  geometry  of  the  combustor  and  re- 
circulation zone  was  the  same  as  in  the  cold  flow  calculations. 

Both  combustion  efficiency  and  wall  static  pressure  distribution 
data  are  available  for  this  configuration,  and  the  results  of 
the  modular  model  calculation  of  these  quantities  are  shown  in 
comparison  with  the  experimental  data  in  Fig.  3.  The  combustion 
efficiency  shown  was  computed  from  the  mass-average  temperature 
at  each  axial  location  in  the  calculation,  using  the  JANNAF 
temperature-rise  combustion  efficiency  definition  [9]  . The 
"ideal"  exit  total  temperature  required  for  this  calculation 
was  obtained  from  the  data  tabulations  provided  to  the  authors 
by  Dr.  R.  R.  Craig  of  AFAPL;  for  these  data,  T™  ( ideal )=2178°K . 

5 

Note  that  in  these  comparisons,  data  points  are  shown  for  the 
fuel/air  ratio  0.053  value  at  which  the  computations  were  carried 
out.  These  data  points  were  obtained  from  plotted  values  of 
combustion  efficiency  vs  fuel-air  ratio  presented  in  Ref.  [8]  . 

The  band  shown  for  each  data  point  represents  the  range  of  observed 
combustion  efficiencies  as  a function  of  fuel-air  ratio  and  is 
a better  indication  of  the  overall  trend  of  the  combustion  efficiency 
vs  fuel-air  ratio  data  than  are  the  individual  data  points 
themselves . 

It  is  evident  from  the  comparison  shown  in  Fig.  3 that 
the  model  overpredicts  the  observed  combustion  efficiency, 
particularly  for  x/D<3.  In  evaluating  this  comparison,  and  the 
others  to  be  described  in  this  section,  it  should  be  noted  that 
the  combustion  efficiency  data  shown  for  X/D=l  and  X/D=2  were 
obtained  in  different  experiments  using  a shortened  combustor. 

Since  the  cold  flow  L/D=3  results  indicate  that  a recirculation 
zone  length  of  1.55  diameters  (obtained  from  Drewry ' s cold  flow 
data  [l0]  for  the  step  height  used  in  the  experiments  of  Craig 
et  al.  [8]  ) is  appropriate  for  a combustor  L/D=3,  it  is  clear 
that  the  flowfield  in  the  L/D=l  and  L/D=2  combustors  is  significantly 
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different  than  that  assumed  in  the  present  analytical  model. 

The  computed  results  shown  in  Fig.  3 show  a discontinuity  in  the 
combustion  efficiency  vs  axial  distance  curve:  this  discontinuity 
occurs  at  the  end  of  the  recirculation  zone  where  the  transport 
of  energy  and  species  abruptly  becomes  zero  (as  appropriate  for 
the  isoenegetic  wall  boundary  condition  assumed  in  the  model). 

The  overprediction  of  the  combustion  efficiency  for  X/D<2 
results  in  significant  differences  in  the  computed  and  measured 
wall  pressure  distributions,  as  is  also  shown  in  Fig.  3.  This 
discrepancy  tends  to  support  the  combustion  efficiency  comparison 
irrespective  of  the  differences  between  the  computational  and 
actual  combustor  configurations  as  outlined  in  the  preceding 
paragraph;  that  is,  the  sharp  pressure  reduction  observed  in 
the  computational  results  icr  X/D<1  is  consistent  with  an  over- 
prediction of  combustion  efficiency  in  this  region.  Here  it 
should  be  noted  that  while  in  the  computational  results  the 
wall  static  pressure  is  normalized  with  respect  to  the  total 
pressure  at  the  dump  plane  the  data  are  normalized  with  respect 
to  a total  pressure  measured  somewhat  upstream  of  the  dump  plane. 

A major  reason  for  the  overprediction  of  combustion  efficiency 
shown  in  Fig.  3 is  the  use  of  a simple,  one-step  global  model 
to  characterize  the  combustion  process.  Since  it  is  known  from 
other  work  that  the  one-step  chemistry  model  used  in  this 
calculation  significantly  underpredicts  the  ignition  delay  time 
and  overpredicts  the  overall  reaction  rate,  further  computations 
were  carried  out  with  modifications  to  the  one-step  global  model 
to  assess  the  effects  of  the  chemistry  model  on  these  reacting 
flow  calculations. 

Results  of  two  additional  computations  of  the  reacting 
sudden  expansion  combustor  flowfield  configuration  are  shown  on 
Fig.  4,  with  the  original  results  (denoted  model  I)  included 
for  comparison.  For  the  calculations  shown  as  model  II,  the 
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pre-exponential  coefficient  in  the  global  model  was  arbitrarily 
reduced  for  X/D<1.55,  returning  to  its  original  value  for 
X/D>1.55;  the  model  used  in  the  well-stirred  reactor  recirculation 
zone  module  was  unchanged.  While  this  modification  has  the  desired 
effect  of  reducing  the  predicted  combustion  efficiency,  and  in- 
creasing the  computed  wall  static  pressure,  the  shapes  of  both 
the  combustion  efficiency  and  wall  static  pressure  curves  are 

X 

incorrect,  particularly  for  g <2.  Thus  an  alternate  approach 
was  taken  for  the  model  III  calculations:  the  original  one-step 
model  was  used  in  the  directed  flow,  while  the  pre-expoential 
coefficient  for  the  recirculation  zone  was  decreased.  As  can 
be  seen  from  Fig.  4,  this  model  produced  a considerable  improve- 
ment in  both  the  shape  and  level  of  the  combustion  efficiency 
curve  for  X/D<2 , as  well  as  in  the  wall  static  pressure  distribution, 
at  the  cost  of  some  accuracy  in  the  prediction  of  combustion 
efficiency  at  X/D=3. 

Because  of  the  fast  kinetics  inherent  in  the  one-step  model 
used  in  the  directed-f low  portion  of  the  model  III  calculations, 
the  flowfield  downstream  of  the  recirculation  zone  (X/D>1.55)  is 
effectively  mixing  limited.  In  all  of  the  computations  described 
to  this  point,  the  two-equation  turbulence  model  has  been  used 
with  coefficients  that  are  essentially  those  recommended  by 
Launder  et  al.  [7]  , the  exception  being  the  exchange  coefficient 
for  the  dissipation  rate,  Pr£ , which  was  taken  to  be  1.0  instead 
of  the  value  of  1.3  used  by  Launder,  et  al.  These  coefficients 
have  been  established  for  use  in  nonreacting  flowfield  calculations, 
and  other  researchers  have  suggested  the  use  of  slightly  different 
coefficients  for  reacting  flowfields.  To  investigate  the  effects 
of  the  turbulence  model  coefficients  on  the  reacting  flowfield 
prediction  using  the  present  model,  computations  were  carried 
out  using  the  coefficients  suggested  by  Hutchinson,  et  al.  [ll]  : 

Cy  = 0.09,  C£l=1.44,  C£2=l . 92 , Prk=0.90,  Pr£=1.22, 
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which  can  be  compared  to  the  coefficients  used  in  the  preceding 
computations : 

0^  = 0. 09,  Ce1=1.40,  Ce2=1.95,  Prk=1.00,  Pre=1.00, 

from  which  it  can  be  seen  that  the  most  substantial  changes  are 
in  the  kinetic  energy  and  dissipation  rate  exchange  coefficients. 
These  differences  will  reduce  the  diffusion  of  dissipation  rate, 
e,  while  increasing  that  of  the  turbulent  kinetic  energy.  As 
the  results  shown  in  Fig.  5 indicate  the  net  effect  of  the 
differences  between  turbulence  model  coefficients  is  to  increase 
the  directed  flow  mixing  rate,  as  is  evidenced  both  by  the  in- 
creased level  of  combustion  efficiency,  and  by  the  change  in 
the  predicted  static  pressure  distribution.  On  the  other  hand, 
the  small  differences  between  the  two  turbulence  model  predictions 
of  combustion  efficiency  in  the  recirculation  zone  region, 

X/D<1.55,  can  be  taken  to  indicate  that  the  changes  inherent  in 
the  use  of  the  turbulence  model  coefficients  of  Hutchinson,  et 
al.  [ 1 l]  have  little  effect  on  the  overall  transport  across  the 
dividing  streamline  shear  layer. 

Also  shown  on  Fig.  5 is  the  effect  of  decreasing  the  effective 
dissipation  length  scale  in  the  recirculation  zone  on  the 
combustion  efficiency  and  wall  static  pressure  distribution 
calculations:  as  in  the  preceding  comparison,  these  computations 

use  the  model  III  global  chemistry  formulation,  with  a reduced 
pre-expoential  rate  coefficient  within  the  sitrred  reactor  model 
of  the  recirculation  zone.  In  the  cold  flow  calculations,  a 
change  of  effective  recirculation  zone  length  scale  l^z  from 
IrZ  = k/2  to  lpZ  = h/4  had  a substantial  effect  on  the  wall 
static  pressure  distribution,  as  seen,  for  example,  in  Fig. 

2f,  with  decrease  in  the  effective  length  scale  resulting  in 
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This  effect  is  still  apparent  in  the  reacting  flow  results  in 
Fig.  5,  but  is  substantially  reduced,  as  is  clear  from  a comparison 
of  the  results  shown  in  Fig.  5 with  those  shown  in  Fig.  2f. 


Although  the  results  shown  in  Figs.  4 and  5 indicate  a 
reasonably  good  level  of  agreement  with  the  data  of  Craig,  et 
al.,  particularly  in  view  of  the  expanded  scale  used  for  the  static 
pressure  comparisons,  all  of  the  predictions  show  a characteristic 
static  pressure  decline  followed  by  a region. of  relatively  uniform 
and  then  slowly  declining  pressure.  This  behavior  does  not  seem 
to  appear  in  the  experimental  data.  However,  other  data  obtained 
by  Craig  using  the  same  combustor  configuration  does  show  the 
characteristic  behavior  just  described.  These  data  were  obtained 
for  premixed  fuel  and  air  in  a combustor  identical  to  that  used 
for  the  comparisons  so  far  described  (L/D=3,  A^/ A^2 .25)  except 
that  an  area  ratio  0.50  convergent  nozzle  was  used  instead  of 
the  A*/Ag=0.40  nozzle  used  to  obtain  the  experimental  results 
previously  shown.  From  the  calculational  standpoint,  the 

differences  between  the  two  sets  of  data  involve  the  inlet  Mach 
number  and  mass  flow  rate.  To  determine  whether  the  typical 
pressure  ratio  behavior  observed  in  the  previous  results  would 
also  be  observed  for  the  initial  conditons  appropriate  to  the 
A*/Ag=0.50  data,  a computation  was  made  using  as  initial  conditions 

an  inlet  total  temperature  of  554°K  (997°R),  an  inlet  static 

2 

pressure  of  1.83  atm  (3881  lb/ft  ),  an  inlet  mass  flow  rate  of 
1.80  kg/sec  (3.45  lb/sec)  and  an  inlet  Mach  number  of  0.401. 

The  model  III  global  kinetics  formulation  was  used,  with  the 
turbulence  model  coefficients  of  Hutchinson,  et  al.  [11]  and 
IrZ  = h/2.  Results  of  this  computation  are  shown  in  Fig.  6. 

In  the  absence  of  data  on  combustion  efficiency  as  a function 
of  L/D  for  this  configuration  the  data  bands  for  the  A*/Ag=0.40 
results  are  again  shown,  for  comparison.  The  important  observation 
is  that  the  characteristic  wall  static  pressure  behavior  observed 
in  the  A*/A2=0.40  comparisons  is  again  observed,  but  for  this 
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A*/Ag=0.50  case  it  is  in  remarkably  good  agreement  with  the 
experimental  data. 


4.3  SUMMARY  OF  PARAMETER  VARIATION  STUDY 


The  results  of  these  calculations  show  that,  while  the 
introduction  of  the  two-equation  turbulence  model  results  in 
two  additional  coefficients  compared  to  the  model  described  in 
Ref.  6 in  which  a one-equation  turbulence  model  was  used,  the 
resulting  predictions  show  a markedly  reduced  sensitivity  to 
the  values  chosen  for  the  dividing  streamline  shear  layer  initial 
thickness  and  spread  rate  coefficients.  Cold  flow  calculations 
show  that  the  most  sensitive  coefficients  are  those  which  describe 
the  shear  stress  level  along  the  dividing  streamline,  CppZ  > anc* 
the  turbulent  intensity  level  and  turbulence  dissipation  length 
scale  within  the  recirculation  zone,  kRZ  and  e^,  respectively. 
However,  although  k^  and  CppZ  have  been  treated  as  independent 
parameters  in  the  work  reported  herein,  they  can  be  related  to 
each  other  based  on  physical  arguments.  This  coupling  will  be 
investigated  in  future  work  on  the  modular  formulation. 


For  reacting  flow,  the  comparisons  show  that,  at  least  for 
the  data  used  for  the  present  work,  the  modeling  of  the  completeness 
of  reaction  within  the  recirculation  zone  region  has  a strong 
effect  on  the  predicted  results.  In  this  work,  using  a very 
simple  one-step  global  finite  rate  model  in  both  the  directed 
flow  and  recirculation  zone  regions,  good  results  were  achieved 
only  when  the  recirculation  zone  reaction  rate  was  reduced  by 
several  orders  of  magnitude.  For  the  cases  in  which  good 
results  were  achieved,  a recirculation  zone  combustion  efficiency 
of  50-60%  is  indicated. 


With  the  chemistry  model  used  in  the  present  computations, 
the  flowfield  downstream  of  the  recirculation  zone  is  essentially 
mixing  limited.  However,  alteration  of  the  turbulence  model 
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coefficients  from  the  original  set,  developed  for  nonreacting 
flowfield  calculations,  to  a set  recommended  for  reacting  flow- 
field  computations,  produced  little  change  in  the  overall  results 
Similarly,  the  sensitivity  of  the  computations  to  the  value 
assumed  for  the  recirculation  zone  dissipation  length  scale  was 
considerably  reduced  in  reacting  flow  compared  to  nonreacting 
flowfield  calculations. 

Overall,  the  results  obtained  with  the  modular  model  in- 
corporating a two-equation  turbulence  model  formulation  are 
encouraging.  Calculations  of  reacting  flow  carried  out  using 
constants  established  in  a geometrically  similar  cold  flow 
configuration  show  a reasonable  level  of  agreement  with 
reacting  flow  data.  The  inaccuracies  observed  in  the  reacting 
flow  prediction  can  be  related  to  either  the  inability  of  a one- 
step  global  chemistry  model  to  accurately  predict  the  overall 
progress  of  the  chemical  reactions  occuring  in  a combustor  flow, 
or  to  the  effects  of  a significant  level  of  unmixedness  in  the 
stirred  reactor  region.  This  role  of  numixedness  in  the  stirred 
reactor  region  will  be  addressed  in  future  work  on  modular  model 
development . 


May  1979 


SAI-79-013-WH 


5 . FUTURE  WORK 

Work  on  the  development  and  exploitation  of  the  modular 
model  will  continue,  with  emphasis  on  the  continuing  development 
of  characterization  of  the  turbulent  mixing  and  chemical  kinetics 
in  sudden  expansion  combustor  flowfields.  It  is  worth  noting 
in  this  regard  that  the  reacting  flow  calculations  presented 
here  (which  use  a global  finite  rate  kinetics  model)  averaged 
approximately  l£  minutes  run  time  on  an  IBM  3033  computer; 
three  to  four  iterations  were  involved  to  obtain  convergence  in 
the  stirred  reactor-directed  flow  iteration.  Thus  the  inclusion 
of  fast  routines  utilizing  more  detailed  kinetics  formulations 
is  of  interest  in  order  to  maintain  the  attractive  run  time 
characteristics  of  the  modular  model. 

Close  cooperation  with  Drs.  F.  D.  Stull  and  R.  R.  Craig 
of  AFAPL  will  continue  to  be  maintained  in  order  to  analyze 
and  interpret  AFAPL  sudden  expansion  combustor  data  using  the 
modular  model.  As  was  noted  earlier  in  this  report,  the  results 
presented  here  make  use  of  unpublished  AFAPL  static  pressure 
distribution  data  made  available  to  the  authors  by  Dr.  Craig; 
these  data  are  invaluable  in  the  modular  model  development. 

In  addition  to  further  development  and  exploitation  of  the 
modular  model,  further  work  in  this  area  will  include  the  develop- 
ment of  a spray  combustion  model,  utilizing  the  two-phase  flow 
equations  for  the  detailed  description  of  the  spray  combustion 
process.  Development  of  a unified  (elliptic)  formulation, 
taking  advantage  of  the  modeling  techniques  which  have  proven 
efficacious  in  the  modular  formulation,  will  also  continue. 
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